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ABSTRACT 

We use a combination of analytical theory, numerical simulation, and data 
analysis to study the propagation of acoustic waves along coronal loops. We 
show that the intensity perturbation of a wave depends on a number of factors, 
including dissipation of the wave energy, pressure and temperature gradients in 
the loop atmosphere, work action between the wave and a flow, and the sensi- 
tivity properties of the observing instrument. In particular, the scale length of 
the intensity perturbation varies directly with the dissipation scale length (i.e., 
damping length) and the scale lengths of pressure, temperature, and velocity. 
We simulate wave propagation in three different equilibrium loop models and 
find that dissipation and pressure and temperature stratification are the most 
important effects in the low corona where the waves are most easily detected. 
Velocity effects are small, and cross-sectional area variations play no direct role 
for lines-of-sight that are normal to the loop axis. The intensity perturbation 
scale lengths in our simulations agree very well with the scale lengths we mea- 
sure in a sample of loops observed by TRACE. The median observed value is 
4.35 X 10^ cm. In some cases the intensity perturbation increases with height, 
which is likely an indication of a temperature inversion in the loop (i.e., tempera- 
ture that decreases with height). Our most important conclusion is that thermal 
conduction, the primary damping mechanism, is accurately described by classical 
transport theory. There is no need to invoke anomalous processes to explain the 
observations. 
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Introduction 



Recent observations of oscillations in coronal loops have inspired a great deal of interest 
in the emerging field of "coronal seismology." See, for example, the proceedings of the SOHO- 
13 Workshop (Erdelyi, Ballester, & Fleck 2004) and the NATO Workshop on Turbulence, 
Waves and Instabilities (Erdelyi et al. 2003). The aim of coronal seismology is to uncover 
the properties of the corona by studying their influence on waves. These properties include 
the temperature and density structure, the magnetic field structure, and the magnitudes of 
important transport coefficients. 

Much of the work so far has been concerned with standing waves — either transverse 
standing waves that are revealed as oscillatory displacements of coronal loops and are believed 
to be excited by large flares (e.g., Schrijver, Aschwanden, & Title 2002) or longitudinal 
standing waves that are revealed as Doppler shifts of very hot (> 6 MK) spectral lines and 
believed to be excited by small flares (e.g., Wang et al. 2002). Both types of standing waves 
are observed to decay on time scales of order 1000 s. Nakariakov et al. (1999) and Ofman & 
Aschwanden (2002) concluded that such rapid damping of transverse oscillations implies a 
dissipation by shear viscosity or electrical resistivity that is many orders of magnitude larger 
than classical transport theory predicts (e.g., Braginskii 1965). In contrast, Ofman & Wang 
(2002) and Mendoza-Briceno, Erdelyi, & Sigalotti (2004) concluded that classical thermal 
conduction is adequate to damp the longitudinal waves. 

The work we present here examines the possibility of using propagating longitudinal 
waves for coronal seismology. Specifically, we consider acoustic waves traveling along coronal 
loops. It is commonly understood that the amplitude of such waves is decreased by dissi- 
pation of the wave energy (damping), but it is not widely appreciated that the amplitude 
can change for a variety of other reasons. Furthermore, these additional effects can actually 
cause the amplitude to increase as the wave propagates. As a simple illustration, consider 
an acoustic wave propagating along a static coronal loop in the absence of dissipation. The 
energy flux of the wave is then conserved, so 
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is the wave's energy density. Here, 5v is the velocity perturbation amphtude, and P, p, and 
T are the pressure, mass density, and temperature, respectively. Making use of the ideal gas 
law, P — RpT, we obtain 

Sv oc -7^, (3) 

{PAf^ 

which shows that the velocity amplitude will vary as the wave propagates whenever T, P, 
or A arc nonuniform along the loop. An upwardly propagating wave typically experiences 
increasing temperatures and decreasing pressures, and both act to increase the velocity 
amplitude. This is not a new result. The direct influence of gravitational stratification 
on acoustic wave propagation has been studied by Ofman, Nakariakov & Deforest (1999), 
Nakariakov et al. (2000), Ofman, Nakariakov, & Sehgal (2000), and De Moortel & Hood 
(2004). Stratification can also have an indirect effect by modifying the dissipation rate 
(Mendoza-Briceno, Erdelyi, & Sigalotti 2004). 

The assumption of no dissipation in the above example is unrealistic. Linear theory 
indicates that acoustic waves propagating in the corona will be significantly damped, pri- 
marily by thermal conduction and to a lesser extent by compressive viscosity and radiation 
(Porter, Klimchuk & Sturrock 1994; Dc Moortel & Hood 2003). Unfortunately, the theory 
has limited validity in the transition region and lower corona, where the temperature and 
pressure can deviate significantly over the distance of a wavelength. To study wave damping 
in this situation requires a full nonlinear treatment (i.e., numerical simulation). We here con- 
sider nonlinear damping by thermal conduction in a non- isothermal loop environment. See 
Ofman, Nakariakov, & Sehgal (2000) for a treatment of nonhnear damping by compressive 
viscosity in an isothermal plume environment. 

Equation 3 concerns the velocity perturbation of the wave, yet it is the intensity pertur- 
bation that is more typically observed. Small- amphtude (4.1 ± 1.5 %) intensity disturbances 
are commonly seen propagating upward along the lower legs of long coronal loops. They 
have been detected by TRACE (Nightingale, Aschwanden, & Hurlburt 1999; Berghmans & 
Clette 1999; De Moortel, Ireland, & Walsh 2000; Robbrecht et al. 2001; De Moortel et al. 
2002) and by the Extreme ultraviolet Imaging Telescope (EIT) on the Solar and Heliospheric 
Observatory (SOHO) (Berghmans & Clette 1999). Similar disturbances have been observed 
in polar plumes (Ofman et al. 1997; Deforest & Gurman 1998). Although the disturbances 
could be transient fiows induced by impulsive heating (Reale et al. 2000), it is far more likely 
that they are traveling acoustic waves. They are periodic over several cycles, with periods of 
282 ± 93 s, and they have propagation speeds of 122 ± 43 km s~^ that are comparable to the 
expected sound speed (De Moortel et al. 2002). Alfven waves can be ruled out because they 
would be much faster and because they are incompressible in the linear regime and would 
not produce the observed intensity perturbations. 
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The goal of our study is to evaluate whether these or other acoustic wave observations 
are useful for coronal seismology. Can we infer the thermodynamic and geometric structure of 
loops by measuring the variation of the intensity perturbations? Can we determine whether 
thermal conduction is anomalously large in the corona, as appears to be the case for shear 
viscosity? The answer has important implications for the microphysics of the Sun and of 
astrophysical plasmas in general. 

Our study involves several interrelated parts, and the paper is organized accordingly. 
We begin in Section 2 with an analytical treatment of acoustic wave propagation in a generic 
loop for which the tcmpcratiuc, pressure, and velocity structure is arbitrary. We derive an 
expression relating the scale length of the intensity perturbation to the scale lengths of the 
physical loop parameters. In Section 3, we apply the results to three specific equilibrium loop 
models and determine which aspects of the models most influence the intensity perturbation. 
We then simulate the propagation of waves in the three loops in Section 4. By combining 
the numerical and analytical results, we are able to deduce the classical thermal conduction 
damping length throughout the loop, including those regions where linear theory breaks 
down. In Section 5, we reanalyze the TRACE observations of De Moortel et al. (2002) and 
compare them directly with the wave simulations. We conclude with some final remarks in 
Section 6. 



2. Analytical Treatment 

As we have indicated, the observed intensity perturbation of an acoustic wave can evolve 
because of dissipation of the wave energy, nonuniformity of the atmosphere through which 
the wave propagates, and work action between the wave and a flow. In addition, the detected 
intensity will depend on the properties of the observing instrument. If coronal seismology is 
to be a useful diagnostic tool, then we must isolate these different influences to determine 
which are the most important. 

For a line-of-sight passing orthogonally through the axis of a coronal loop, the observed 
intensity of the optically-thin plasma is 

I{s) = n^G{T)A^'\ (4) 

where n is the electron number density, A is the cross-sectional area, G(T) is the instrument 
response function describing the temperature-dependent sensitivity of the observing instru- 
ment, and s measures position along the loop. It is assumed here that the cross section 
is circular, as suggested by both observational (Klimchuk 2000) and theoretical (Khmchuk, 
Antiochos, & Norton 2000) considerations, and that the plasma is uniform over the cross 
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section. See Ofman, Nakariakov, & Sehgal (2000) for a discussion of the effects of cross-field 
gradients. 

An acoustic wave propagating along the loop will perturb the density and temperature 
and therefore the intensity. The cross section will be nearly unaffected, however, because 
/3 = SttP/B'^ is very small and the magnetic field behaves hke a rigid tube. Differentiating 
Equation 4 gives 

51 = 2nG(T)A^/^5n + n''A^/^^5T, (5) 

where 51, 5n, and 5T are the intensity, density, and temperature perturbation amplitudes, 
respectively. As discussed in the Appendix, the second term on the right-hand-side is ex- 
pected to be considerably smaller than the first term for waves observed by TRACE. We 
therefore drop the second term to obtain the approximate expression 

51 ^ 2nG{T)A^/^5n. (6) 

It is interesting that the temperature perturbation has no direct effect on 51 even though it 
is crucially important for thermal conduction damping. 

In the absence of dissipation, the energy density of an acoustic wave satisfies 

C 

e oc (7) 

where v is the velocity of a flow that may be present in the loop (e.g., Jacques 1977, Equation 
19). This form accounts for the exchange of work energy between the wave and the flow. It 
reduces to the familiar expression given in Equation 1 when v = 0. We can add the effect of 
dissipation by multiplying the right-hand-side of Equation 7 by e~^^^'^, where 

is the dissipation scale length for the wave energy. A negative sign is used so that Ha will be 
positive for waves propagating in the direction of increasing s, to follow usual convention. We 
refer to as the "damping length" and stress that it refers exclusively to energy dissipation 
by direct plasma heating. Energy loss by work is not included in Ha- 

Comparing the modified form of Equation 7 with Equation 2, we find that the velocity 
perturbation amplitude is given by 

1 rc\^^ 



5v oc {^] e ^"d. (9) 
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This is in agreement with the findings of De Moortel & Hood (2004) for the case of uniform 
temperature and no flow. From the continuity equation, the velocity and density perturba- 
tions are related according to 

Combining Equations 6, 9, and 10, we obtain 

{Cs + v)Cl'^ 

Noting that Cg oc T^l"^ and n oc P/T, and deflning the scale length of any quantity / to 
be Hf = — f / {df / ds) , we finally arrive at an expression for the scale length of the intensity 
perturbation: 



Hsr^ 

where M = v/Cg is the Mach number, and 



(12) 



Tg ^ G{T) (^^) (13) 

is the scale temperature of the instrument response function. Equation 12 is very useful 
because it reveals how the different properties of the loop affect the intensity perturbation. 
It shows that the characteristic distance over which 51 varies depends directly on the char- 
acteristic distances over which the loop parameters vary. The dependence is linear in the 
sense that H^i is proportional to the parameter scale length when only one parameter is 
non-constant. For example, H^j = [2/3)Hp with pressure stratification alone. 

The terms with the smallest scale lengths in Equation 12 have the biggest affect on Hsi- 
Furthermore, the smaller a scale length is, the smaller it makes Hsi- Consider the case of 
a static loop {v = 0). If we treat each term separately, then H^j oc {2/3)Hp, {Tg/T)Ht, 
{—4:/9)Ht, or 2Hd. We indicate two separate temperature terms because one depends on 
the observing instrument and the other does not. The coefficients of the terms are all of 
order unity^ , so the relative importance of the terms is determined primarily by the relative 
sizes of the scale lengths. 

The presence of a flow generally has minimal impact on Hgi. The (1 -|- M)H^/M term 
tends to be large (i.e., not important) because velocities tend to be small where velocity gra- 
dients are large, and vice versa. The change to the coefficient of the instrument-independent 
Ht term is minimal, with its magnitude increasing only from 4/9 to 1/2 for M — 1. 



^Ta/T 1/3 for spectral line observations and can be much larger for broad-band observations. 
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Since the terms in Equation 12 can be either positive or negative, the intensity pertur- 
bation can either decrease or increase as the wave propagates. Upwardly propagating waves 
generally experience decreasing pressures and increasing temperatures, and such conditions 
cause the intensity perturbation to decrease with height. This may seem inconsistent with 
our demonstration in the Introduction that the velocity perturbation increases with height 
under normal conditions. In fact, there is no contradiction. Intensity and velocity pertur- 
bations are not expected to behave in the same way, which is a subtle but important point 
that is sometimes overlooked. Prom Equations 3, 6, and 10, we have that 5v oc rV^p-i/s 
but that 51 oc j'-9/4p3/2^ ^j^g velocity and intensity perturbations depend on temperature 
and pressure in an opposite sense. 

It is interesting that H^i has no direct dependence on the cross-sectional area of the loop. 
This is only true when the line-of-sight is orthogonal to the loop axis, as we have assumed, 
because then the area dependence of Equation 4 exactly cancels the area dependence of 
Equation 7. There is nonetheless an indirect dependence on the cross section due to its 
effect on the temperature and pressure structure of the loop atmosphere. The velocity 
perturbation, unlike the intensity perturbation, does depend directly on the cross section, as 
shown in Equation 9 and demonstrated in De Moortel & Hood (2004). 

In summary, the evolution of the intensity perturbation of a propagating acoustic wave 
depends primarily on the pressure and temperature structure of the loop and the wave energy 
damping length. The relative importance of these different factors depends on the details of 
the loop atmosphere. In the next section, we evaluate the separate terms in Equation 12 for 
three specific equilibrium loop models. 



3. Equilibrium Loop Models 

In the solar corona, where the electrical conductivity is very large and j3 is very small, the 
plasma is confined by the magnetic field and its evolution is described by the one-dimensional 
conservation equations of mass, momentum and energy: 

''"'^A'-^-'Z-m-^, (15) 



dt A ds ds 

dE ld[A{E + P)v] ^df^ rr.5/29T\ 2a/^N ^ 

-m + A OS = + AO's (^'^"^ ' ) - + ^' 



(16) 



where 

1 9 ^ 
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As before, T, P, p, v, and A are the temperature, pressure, mass density, velocity, and 
cross-sectional area, respectively; s is distance along the loop from the left base of the 
model; g\\{s) is the component of gravity parallel to the loop axis; Q{s) is the volumetric 
heating rate; A(T) is the optically-thin radiation loss function; k,q = 10~^ is the coefficient 
of thermal conduction; and 7 = 5/3 is the ratio of the specific heats. We assume a fully 
ionized hydrogen plasma, so p = rripn and P = 2kBnT, where n is the electron number 
density and is Boltzmann's constant. We do not include compressive viscosity, because 
it is generally small under the conditions we will investigate (Porter, Klimchuk & Sturrock 
1994; De Moortel & Hood 2003). Finally, a — a(s, t) is a spatially locahzed acceleration that 
we impose to generate acoustic waves. It is set to zero for the equilibrium solutions. 

The geometry of the loop is defined by the specific form of g\\ (s) . We adopt a semicircular 
shape for the coronal part, so that 

2 

cos [0{s)] , (18) 

where h{s) is the height above the solar surface, Qq is the gravitational acceleration at the 
surface, Rq is the solar radius, and d{s) is the angular position along the semicircle. This 
form accounts for the decrease in gravity with height, which can be important for high- 
arching loops. We choose a coronal length L = 3 x 10^'' cm. Attached to the corona at 
both ends are 6 x 10^ cm long chromospheric sections, giving a total length for the model of 
4.2 X 10^° cm. 

The optically thin radiative loss function, A(T), is a piecewise continuous power law as 
described in Klimchuk & Cargill (2001). It is based on the atomic physics calculations of J. 
Raymond (private communication) and uses abundances twice as large as the coronal values 
given by Meyer (1985). The loss function drops precipitously to zero between 3.0 x 10^ and 
2.95 X 10^ K in order to maintain a nearly constant chromospheric temperature within this 
range. The detailed physics of the chromosphere involves optically-thick radiative transfer 
and is not treated, but our model accurately reproduces the energetic and dynamic coupling 
between the corona and chromosphere which is vital to a realistic simulation. 

Although there is evidence that many coronal loops are inherently time dependent (e.g., 
Cargill & Klimchuk 1997), we limit ourselves in this study to equilibrium loops. We consider 
two static models and one model with steady end-to-end fiow. Static equilibrium requires 
that the heating function, Q{s), be symmetric with respect to the loop midpoint. In order to 
produce as flat a temperature profile as possible, as suggested by Tii^Ci? observations (Lenz 
et al. 1999; Aschwanden, Nightingale, & Alexander 2001), we adopt a heating function 



9\\{s)^gQ 



Rg) 



Rq + h(s) 
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that decreases exponentially with height in both legs. In the "left" leg, 

Q{s) = Qo exp (-^) , (19) 

where Hq = 7.5 x 10^ cm is the heating scale length and Qq = 4.6 x 10~^ erg cm~^ s^^ is the 
heating rate at the left footpoint (sq = 6 x 10^ cm). The heating function is a mirror image 
in the "right" leg. We have set the magnitude of the heating to produce a peak temperature 
of approximately 1.2 MK, characteristic of Ti^^lC^^ loops. 

To obtain a steady flow solution, we follow Patsourakos, Klimchuk & MacNeice (2003) 
and apply an exponential heating function over the entire loop, without mirror reflection. 
The heating rate decreases from the left base, past the loop midpoint, all the way to the 
right base. In this case we use a scalelength Hq — L/10 — 3 x 10^ cm and magnitude 
Qo = 6.5 X 10~^ erg cm~^ s~^, again chosen to produce a peak temperature near 1.2 MK. 
The extreme heating asymmetry drives a strong left-to-right flow in the loop. 

Coronal loops that are visible along their entire length tend to have constant cross 
sections (Klimchuk et al. 1992; Klimchuk 2000; Watko & Klimchuk 2000; Lopez-Fuentes & 
Klimchuk 2004). Acoustic waves are typically observed in much longer loops, for which only 
the lower legs are visible. There is some evidence that these loops expand with height. De 
Moortel et al. (2002) measured an average expansion rate of 0.28 ± 0.16 (width change per 
unit length). The average width of the loops is 8.1 ± 2.8 x 10^ cm, so the expansion rate 
can be expressed as 3.5% per 10^ cm. In one of our models, we assume a cross section that 
expands according to 



A{s) = A^id COS 



(20) 



2 ^mid 

where Amid is the area at the loop midpoint {smid = 2.1 x 10^° cm) and we set x — 1-3 to give 
a width expansion rate of 4.3% per 10^ cm at the base of the corona (area expansion rate 
of 8.7%). The cross section is taken to be constant throughout most of the chromosphere. 
Figure 1 shows the variation of area with position in the loop. 

Equations (14-17) are solved numerically using a code we call Adaptively Reflncd Go- 
dunov Solver (ARGOS). As described in Antiochos et al. (1999), it utilizes a second order 
Riemann type solver (Van Leer 1979) to treat the one-dimensional equations in the absence 
of sources. The sources arc accounted for using an operator splitting method in second 
order. An important feature of ARGOS is the parallel adaptive mesh refinement package 
PARAMESH, which has been set to refine or derefine the local grid based upon density 
variations. Rigid wall boundary conditions are imposed at the ends of the model, which are 
many gravitional scale heights deep in the chromosphere. To obtain an equilibrium solution, 
we begin with a reasonable guess of the loop atmosphere based on well-known scaling laws 
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(e.g., Rosner, Tucker, & Vaiana 1978), and allow the solution to relax asymptotically to a 
steady state. 

We consider three different equilibrium loops: Case 1 is a constant cross section loop 
with symmetric heating (static equilibrium); Case 2 is a constant cross section loop with 
asymmetric heating (steady flow equilibrium); and Case 3 an expanding cross section loop 
with symmetric heating (static equilibrium). Figures 2 and 3 show the profiles of temperature 
and pressure, respectively, for all three cases. Figure 4 shows the profile of velocity for Case 
2. The maximum velocity is 46 km s^^ and corresponds to a Mach number of 0.29. There 
are clear differences among the temperature and pressure profiles for the three cases. In 
particular, the temperature profile in Case 2 has an inversion that is not present in Cases 
1 and 3. The physical explanation for this behavior is given in Patsourakos, Klimchuk & 
MacNeice (2003). Differences between Cases 1 and 3 can be understood in terms of changes 
to the energy balance associated with the constriction of the loop footpoints. For Case 3, a 
greater fraction of the coronal heating energy must be radiated from the upper parts of the 
loop and less conducted down to the lower parts, so the density is higher and the temperature 
profile flatter. 

Figures 5, 6, and 7 indicate the values of the terms in Equation 12 as a function of 
position along loop for Cases 1, 2, and 3, respectively. Only the left half of the loop is 
shown. The iJp, Ht-, and Hy terms were obtained directly from the equilibrium models. Of 
these, the pressure term is most important (i.e., has the smallest magnitude) throughout 
most of the corona. Pressure is gravitationally stratified and therefore Hp increases toward 
the apex as the loop bends to become more and more horizontal. The temperature terms 
are important in the lower corona, and they dominate in the transition region where the 
gradients are very steep. The interesting shape of the instrument-dependent Ht curve is due 
to the form of G{T). We have used a G{T) corresponding to the TRACE 171 channel, which 
peaks at T = 0.95 MK. The curve is vertical at the position where this temperature occurs, 
because Tg oc {dG/dT)~^ switches discontinuously from large positive to large negative 
values. The Ht terms are negative in the plotted portion of Case 2. This is because of the 
temperature inversion. The H^ term is insignificant. 

4. Wave Simulations 

In principle, we can use linear theory to estimate the value of the damping length from 
the equilibrium models and then solve for H^i in Equation 12. Following (Porter, Klimchuk 
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fe Sturrock 1994, Eqn. 4.1.2f), we have that 



^ 146—, 



(21) 



where the asterisk indicates that the expression apphes to adiabatic waves propagating in a 
uniform medium. For Case 1, fa 9.7 x 10^ cm at position s = 1.0 x 10^'^ cm. We show 
in the Appendix, however, that waves of the observed periods are approximately isothermal 
under the conditions in the loop. Furthermore, the wavelengths are not small compared to 
both the pressure and temperature scale lengths. Equation 21 is therefore not vahd. To 
determine Hsi and we must perform a nonlinear simulation of wave propagation. 

We generate waves in our loops using a sinusoidal acceleration that is localized just 
above the left footpoint: 



The Gaussian spatial profile has a 1/e halfwidth a — 10^ cm and is centered at Si = 6.2 x 10^ 
cm for Cases 1 and 3, and si = 8.5 x 10^ cm for Case 2. The temperature at these locations is 
0.55, 1.20, and 0.73 MK for Cases 1, 2, and 3, respectively. We choose an oscillation period 
T = 250 s, similar to those observed by TRACE, and an acceleration magnitude oq that 
produces a wave amplitude ~ 4% {SI ~ 8%), about twice as large as the average observed 



The left panel of Figure 8 shows the velocity perturbation in Case 1 at a time when six 
wave fronts have been generated. The right panel shows the corresponding intensity pertur- 
bation computed from Equation 4 for a simulated TRACE 171 observation. Its amplitude 
decreases rapidly as the wave propagates along the loop. By following one of the wave fronts, 
we determine the scale length 



The solid curve in Figure 9 shows Hsi as a function of position along the loop for the second 
wave front. The other wave fronts give similar scale lengths. The dashed curve in the figure 
shows Hsi for the 195 channel. Differences are due entirely to differences in temperature 

responses of the two channels. The wiggles in the curves are not physical, but arise from 
small inaccuracies in determining the position of the wave front (the position of maximum 
perturbation) on the finite numerical grid. A smoothed version of the 171 curve is shown as 
the solid curve in Figure 5. Corresponding curves for Cases 2 and 3 are shown Figures 6 and 
7. We see that Hsi is typically between 1 and 6x10^ cm and somewhat larger near the top 




(22) 



by TRACE. 




(23) 
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of the loop in Case 1. This agrees very well with the actual values observed by TRACE, as 
discussed in the next section. 

The dotted curves in Figures 5, 6, and 7 indicate 211^ as computed from Equation 12 
using Hgi from the simulations. The damping term is of comparable importance to the 
pressure and temperature terms in the low corona. It is of significantly greater importance 
in the Tipper corona, although waves there are difficult to detect because of their reduced 
intensity (Fig. 8). We note that is larger than H^. This is expected since (approximately) 
isothermal waves dissipate more slowly than adiabatic waves (De Moortel & Hood 2003). 



Using TRACE 171 observations, De Moortel et al. (2002) studied waves propagating 
upward in the lower corona of long loops. They performed a wavelet analysis on the intensity 
perturbations and measured a quantity they call the "detection length." It is defined to be 
the distance over which the wavelet power exceeds the 99% confidence level with respect to 
photon counting noise (i.e., the distance over which the intensity perturbation is at least a 
2.1(7 detection). At greater altitudes the oscillation signal vanishes into the noise. 

The relationship between the detection length and the intensity perturbation scale 
length is not straightforward. Whereas Hsi indicates the characteristic distance over which 
the intensity perturbation varies and is independent of the signal strength, the detection 
length depends on both the perturbation amplitude and the intensity of the unperturbed 
loop. For a given Hgj, a large amplitude wave would be "detectable" out to greater distances 
than a small amplitude wave and would therefore have a longer detection length. 

To facilitate a direct comparison with our simulations, we measure an observed intensity 
scale length: 



where 5Ii and 5I2 are the intensity perturbation amplitudes at positions si and S2, which are 
the first and last positions where the wavelet power meets the 99% criterian. We estimate 
5Ii and 5I2 as follows. We first compute a time sequence of the intensity at each position 
by averaging over a cross-sectional slice that is 2-4 pixels thick. We then subtract a linear 
fit of the time sequence, i.e., we remove the unperturbed signal taking into account any 
fading or brightening trend that may be present. Finally, we compute a root-mean-square 
(RMS) intensity from the residuals, which we claim is directly proportional to the intensity 
perturbation amplitude. The constant of proportionality is unimportant, since it falls out of 
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Observations 
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the ratio in Equation 24. 

Table 1 lists 6Ii, 6I2, As = S2 — Si (the detection length), and Hgj'^ for the 17 loops 
studied originally by De Moortel et al. (2002). Several of the loops were observed multiple 
times. The distribution of Hgj^ is shown as a histogram in Figure 10. Values greater than 
8 X 10^ cm have been assigned to the 8x10^ cm bin. Nearly one-third of the scale lengths are 
negative. We choose not to plot them to the left of the positive values since a more sensible 
physical ordering is to begin with wave amplitudes that decrease rapidly with height (small 
positive H'^Y), progress to wave amplitudes that decreases slowly with height (large positive 
Hfj'), and finally to wave amplitudes that increase with height (negative H'^Y). We therefore 
lump the negative values together and assign them to the 9 x 10^ cm bin. With these changes, 
the median (not the mean) of the distribution is 4.35 x 10^ cm. The H^i from our wave 
simulations generally range between 1 and 6x10^ cm, which is in excellent agreement with 
the observations. 

The negative Hf}^ measurements, indicating intensity amplitudes that increase with 
height, are most likely a result of a temperature inversion. The pressure term in Equation 
12 is positive due gravitational stratification; the dissipation term is positive by definition; 
and the fiow term is probably very small, as we have discussed. Only the temperature terms 
arc left to make Hf^^ negative. Since < for T > 0.95 MK, Ht must be positive 
and temperature must decrease with height (i.e., there must be a temperature inversion 
compared to the usual situation). We have seen that a temperature inversion is present 
in the steady flow equilibrium of Case 2. The temperature gradient is not steep enough, 
however, to drive H^i negative in this model. De Moortel, Parnell, & Hood (2003) found 
evidence for temperature inversions in two loops of their sample for which they had both 
171 and 195 observations and could apply filter ratio diagnostics. We plan to examine other 
cases for which data in both channels are available. 

Some readers may find it strange that H^Y could be negative at the same time that the 
detection length is finite. This situation occurs whenever the unperturbed loop intensity, and 
hence photon counting noise and detection threshold, increase with height faster than the 
intensity perturbation. The effect of increasing detection threshold with height also explains 
why positive Hg^ values tend to be larger than or much larger than the corresponding 
detection lengths. 

We caution that the uncertainties in our measurements are quite large. Following Klim- 
chuk & Gary (1995), it can be shown that the uncertainty in Hg^ is given approximately 
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where ASIi/SIi and ASI2/5I2 are the fractional uncertainties in the intensity amphtudes. 
It is difficult to determine the formal uncertainties in the intensity amplitudes, but visual 
inspection of the time sequences suggests that 15% is a reasonable estimate. With this 
value, the uncertainties in the scale lengths are as listed in the final column of Table 1. It 
is evident that many of the scale length measurements arc, by themselves, rather dubious. 
However, because the measurement errors are of a random nature, so that the true scale 
length is sometimes underestimated and other times overestimated, and because the sample 
of measurements is reasonably large, we feel that the basic results discussed above are still 
valid. 



6. Discussion 

We have shown that the intensity perturbation of an acoustic wave propagating along 
a coronal loop depends on a number of factors, including dissipation of the wave energy, 
pressure and temperature structure in the loop atmosphere, work action between the wave 
and a flow, and the sensitivity properties of the observing instrument. For equilibrium loops, 
dissipation and pressure and temperature stratification dominate in the low corona where 
the waves are most easily detected. Unfortunately, the effects are of comparable magnitude, 
and it would be difficult to disentangle them to infer the structure of an unknown loop 
atmosphere. Coronal seismology offers no big advantage over spectroscopy in this regard. 

We have also shown that the intensity perturbations observed in 1-2 MK loops by 
TRACE are easily reproduced by propagating wave simulations with classical thermal con- 
duction damping. There is no need to invoke anomalously large conduction or a different 
damping mechanism altogether. Ofman & Wang (2002) and Mendoza-Briceno, Erdelyi, & 
Sigalotti (2004) similarly concluded that classical thermal conduction is adequate to explain 
the damping of longitudinal standing waves observed as Dopplcr shifts in much hotter (> 6 
MK) loops. The "loops" in their simulations are non-radiating isothermal structures that 
are disconnected from the chromosphere and therefore do not include the complex energy 
balance of our simulations, but the basic conclusion concerning the effectiveness of classical 
thermal conduction damping should be correct. 

Thermal conduction should not affect the transverse standing waves that are revealed as 
oscillating loops in Ti^^lC^^ movies (Schrijver, Aschwanden, & Title 2002; Aschwanden et al. 
2002). Nakariakov et al. (1999) and Ofman & Aschwanden (2002) suggested that these waves 
are instead damped by shear viscosity and/or electrical resistivity. However, the magnitude 
of the damping must be at least 10^ times larger than predicted by classical transport theory. 
Our study rules out the possibility that either thermal conduction or compressive viscosity 
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is enhanced by a similar factor. 

It is not unreasonable that a mechanism of anomalous transport would affect the various 
transport processes differently. This is known to be the case for classical transport whenever 
a strong magnetic field is present. Cross- field shear viscosity is greatly inhibited compared 
to field-ahgned compressive viscosity and thermal conduction. The difference in the viscous 
coefficients is 10-12 orders of magnitude under coronal conditions. Anomalous transport 
might also be highly anisotropic, though not necessarily in the same sense. For example, 
turbulence would enhance shear viscosity much more than either compressive viscosity or 
thermal conduction (Braginskii 1965, p. 236). Wc conclude that the damping of transverse 
waves by anomalous shear viscosity and the damping of longitudinal waves by classical ther- 
mal conduction are not inconsistent. However, further detailed study of particle transport is 
called for, since there are many important applications throughout astrophysics and space 
physics. 

We close by pointing out that the equilibrium models used in our study are not neces- 
sarily a good representation of the loops in which propagating waves are observed. Those 
loops are very long and tend to occur at the perimeters of active regions. A majority of more 
typical TRACE loops are known to be inconsistent with both static equihbrium (Aschwan- 

den, Schrijver, & Alexander 2001; Winebarger, Warren, & Mariska 2003) and steady flow 
equilibrium (Patsourakos, Klimchuk & MacNeice 2003), even if the heating is concentrated 
near the footpoints. These loops are best explained as bundles unresolved strands that are 
heated impulsively and quasi-randomly by nanoflares (Cargill 1994; Klimchuk 2002; Warren, 
Winebarger, & Mariska 2003). In this picture, a given loop contains a wide distribution of 
strand temperatures at all times. 

Whether the long loops with waves are also multi stranded is not known. Marsh et al. 
(2003) recently studied one of the cases in the De Moortel et al. (2002) sample and found that 
the intensity oscillations visible in TRACE 171 are simultaneously visible in lines of He I, 
O V, and Mg IX observed by the CDS instrument on SOHO. This suggests that the loop may 
indeed be multi stranded. If so, then waves launched at the base would propagate at different 
speeds in the different strands. One consequence is that the waves would experience damping 
by phase mixing (Heyvaerts & Priest 1983; Ofman & Aschwanden 2002). Another is that 
the wave front would become progressively distorted and appear to broaden with height, 
giving the impression of more damping than is actually present. We plan to examine this 
latter possibility with additional modeling and a more detailed examination of the TRACE 
data. 

We thank Alan Hood for suggesting this problem to us and Spiro Antiochos for helpful 
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discussions. We also thank the referee for useful comments. This work was supported by 
NASA and the Office of Naval Research. 



A. Isothermal versus Adiabatic Waves 



Acoustic waves can behave isothermally, adiabatically, or somewhere in between de- 
pending on their wavelength and the temperature and density of the plasma in which they 
propagate. For isothermal behavior, 5T — 0, and for adiabatic behavior, 

uTl 

5T^(-f-l)T—. (Al) 
n 

The solid line in Figure 11 shows the ratio of the left to right sides of Equation Al as a 
function of position along the loop for the waves simulated in Case 1, described in Section 
4. The ratio is < 1, but not <^ 1, indicating that the waves are roughly isothermal. This 
is consistent with the result of Porter, Klimchuk & Sturrock (1994, Eqn. C6a), which can 
be modified to state that the transition between isothermal and adiabatic behavior occurs 
when 

2^3/2 

r ^ Tcrit = 65 , (A2) 

n 

with shorter periods corresponding to isothermal behavior. The 250s period of the simulated 
waves is about 3 times shorter than Tcru at s = 10^° cm. 

The dashed line in Figure 11 shows the ratio of the second to first terms on the right- 
hand-side of Equation 5. This ratio is of order 10 - 20 % for most of the loop, justifying our 
neglect of the second term in arriving at Equation 6. 



REFERENCES 

Antiochos, S. K., MacNeice, P. J., Spicer, D. S., & Klimchuk, J. A. 1999, ApJ, 512 985 
Aschwanden, M. J., De Pontieu, B., Schrijver, C, J., & Title, A. M., ApJ,206, 99 
Aschwanden, M. J., Nightingale, R. W., & Alexander, D. 2000, ApJ,541, 1059 
Aschwanden, M. J., Schrijver, C, J., & Alexander, D. 2001, ApJ,550, 1036 
Berghmans, D. & Clette, F. 1999, Sol. Phys., 186, 207 
Braginskii, S. I. 1965, Rev. Plasma Phys., 1, 205 



-17- 



CargiU, P. J. 1994, ApJ, 422, 381 

Cargill, P. J., & Klimchuk, J. A. 1997, ApJ, 478, 799 

Dahlburg, R. B., Klimchuk, J. A., & Antiochos, S. K. 2004, ApJ, submitted 
Deforest, C. E. & Gurman, J. B. 1998, ApJ, 501 L217 
De Moortel, I., & Hood, A. W. 2003, A&A, 408, 755 
De Moortel, I., & Hood, A. W. 2004, A&A, 415, 705 

De Moortel, I., Hood, A. W., Ireland, J., & Walsh, R. W. 2002a, Sol. Phys., 209 61 

De Moortel, I., Ireland, J., & Walsh, R. W. 2000, A&A, 355 L23 

De Moortel, I., Parncll, C. E., & Hood, A. W. 2003, Sol. Phys., 215 69 

Erdelyi, R., Ballester, J. L., & Fleck, B. (ed.) 2003, ESA SP-547, Proceedings of the SOHO- 
13 Workshop on Waves, Oscillations, and Small-Scale Events in the Solar Atmosphere 
(Noordwijk: ESA) 

Erdelyi, R., Petrovay, K., Roberts, B., & Aschwanden, M. (ed.) 2003, NATO Science Ser. 
124, Turbulence, Waves and Instabilities in the Solar Plasma (Dordrecht: Kluwer) 

Heyvaerts, J., & Priest, E. R. 1983, A&A, 117, 220 

Jacques, S. A. 1977, ApJ, 215, 942 

Klimchuk, J. A. 2000, Sol. Phys., 193 53 

Klimchuk, J. A. 2002, Inst, for Theoretical Physics Conference on Solar Magnetism and 
Related Astrophysics, ed. G. Fisher & D. Longcope (Santa Barbara: Univ. California) 
(http:/ /onhne.kitp. ucsb.edu/onhne/solar_c02/klimchuk/) 

Klimchuk, J. A., Antiochos, S. K., & Norton, D. 2000, ApJ, 542, 504 

Klimchuk, J. A. & Cargill, P. J. 2001, ApJ, 553, 440 

Klimchuk, J. A. & Gary, D. E. 1995, ApJ, 448, 925 

Klimchuk, J. A., Lemen, J. R., Feldman, U., Tsuneta, S., & Uchida, Y. 1992, PASJ, 44, 
L181 

Lenz, D. D., DeLuca, E. E., Golub, L., Rosner, R., & Bookbinder, J. A. 1999, ApJ, 517, 
L155 



Lopez-Fuentes, M., $ Klimchuk, J. A. 2004, ApJ, in preparation 
Marsh, M. S., Walsh, R. W., De Moortel, L, & Ireland, J. 2003, A&A, 404, L37 
Mendoza-Briceno, C. A., Erdelyi, R., & Sigalotti, L. G. 2004, ApJ, 605, 493 
Meyer, J.-R 1985, ApJS, 57, 173 

Nakariakov, V. M., Ofman, L., DeLuca, E. E., Roberts, B., & Davila, J. M. 1999, Science, 
285, 862 

Nakariakov, V. M., Verwichte, E., Berghmans, D., & Robbrecht, E. 2000, A&A, 362 1151 

Nightingale, R. W., Aschwanden, M. J., & Hurlburt, N. E. 1999, Sol. Phys., 190, 249 

Ofman, L., & Aschwanden, M. J. 2002, ApJ, 576, L153 

Ofman, L., Nakariakov, V. M., & Sehgal, N. 2000, ApJ, 533, 1071 

Ofman, L., Romoh, M., Poletto, G., Noci, C., & Kohl, J. L. 1997, ApJ, 491, Llll 

Ofman, L., Nakariakov, V. M., & Deforest, C. E. 1999, ApJ, 514, 441 

Ofman, L. & Wang, T. 2002, ApJ, 580, L85 

Patsourakos, S., Klimchuk, J. A., & MacNeice, P. J. 2003, ApJ, 603, 322 

Porter, L. J., Klimchuk, J. A., & Sturrock, P. A. 1994, ApJ, 435 482 

Reale, P., Peres, G., Serio, S., DeLuca, E. E., & Golub, L. 2000, ApJ, 535 412 

Robbrecht, E., Verwichte, E., Berghmans, D., Hochedez, J. P., Poedts, S., & Nakariakov, V. 
M. 2001, A&A, 370, 591 

Rosner, R., Tucker, W. H., & Vaiana, G. S. 1978, ApJ, 220, 643 

Schrijver, C. J., Aschwanden, M. J., & Title, A. M. 2002, Sol. Phys., 206, 69 

Van Leer, B. 1979, J. Comput. Phys., 32, 101 

Wang, T. J., Solanki, S. K., Curdt, W., Innes, D. E., & Dammasch, I. E. 2002, ApJ, 574, 
LlOl 

Warren, H. P., Wmcbargcr, A. R., & Mariska, J. T. 2003, ApJ, 593, 1174 
Watko, J. A. & Klimchuk J. A. 2000, Sol. Phys., 193 77 



- 19- 

Winebarger, A. R., Warren, H. P., & Mariska, J. T. 2003, ApJ, 587, 439 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-20- 




10 20 30 40 

Position (10^ cm) 



Fig. 1. — Cross-sectional area versus position along the loop in Case 3. 
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Fig. 2. — Temperature versus position along the loop in Cases 1 (solid), 2 (dashed), and 3 
(dotted). 
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Fig. 3. — Pressure versus position along the loop in Cases 1 (solid), 2 (dashed), and 3 
(dotted). 
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Fig. 4. — Velocity versus position along the loop in Case 2. 
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Fig. 5. — Values of the terms in Equation 12 versus position along the loop for Case 1: Hp 
and Ht terms are from the equilibrium model; Hsi is from the wave simulations; and 2H(i 
is from the equation. 
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Fig. 6. — Values of the terms in Equation 12 versus position along the loop for Case 2: Hp, 
Ht, and Hy terms are from the equilibrium model; Hsi is from the wave simulations; and 
2Ha is from the equation. 
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Fig. 7. — Values of the terms in Equation 12 versus position along the loop for Case 3: Hp 
and Ht terms are from the equilibrium model; Hsi is from the wave simulations; and 2H(i 
is from the equation. 




Fig. 8. — Velocity perturbation (left) and intensity perturbation (right) versus position along 
the loop in Case 1. Intensities were computed assuming observations made in the TRACE 
171 channel. 
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Fig. 9. — Intensity scalelength versus position along the loop in Case 1 for simulated obser- 
vations in the 171 (sohd) and 195 (dotted) channels of TRACE. Wiggles in the curves are 
numerical artifacts. 
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Fig. 10. — Histogram distribution of observed intensity scale lengths. Values greater than 
8 X 10^ cm have been assigned a value of 8 x 10^ cm, and negative values have been assigned 
a value of 9 x 10^ cm. 
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Fig. 11. — Ratio of the left to right sides of Equation Al (sohd) and ratio of the second to 
first terms on the right-hand-side of Equation 5 (dot-dash) versus position along the loop in 
Case 1. 



Table 1. Observed Intensity Scale Lengths 
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Table 1 — Continued 
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Note. — Intensities are in units of counts 
s~^ pix~^, and lengths are in units of 10^ cm. 



